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Finite-volume lattice QCD calculations offer the possibility of extracting resonance parameters 
from the energy-dependent elastic phase-shift computed using the Liischer technique. In this letter, 
as a trial of the method, we report on the extraction of the non-resonant phase-shift for S and 
Z)-wave 7T7T isospin-2 scattering from dynamical lattice QCD computations. We define a variational 
basis of operators resembling pairs of pions of definite relative momentum and extract a spectrum of 
excited states that maps to phase-shifts at a set of discrete scattering momenta. Computations are 
performed with pion masses between 400 and 520 MeV on multiple spatial volumes. We observe no 
significant quark mass dependence in the phase-shifts extracted which are in reasonable agreement 
with the available experimental data at low momentum. 



Introduction: The hadron spectrum and interactions 
of QCD can be studied from first principles using nu- 
merical simulation of the quark and gluon fields on a fi- 
nite lattice. While significant progress has been made in 
studying isolated excited meson states with qq-like op- 
erators [TJ |5] , it remains challenging to extract proper- 
ties of resonances that appear in the scattering of stable 
hadrons. One procedure, due to Liischer [3], maps the 
discrete spectrum of eigenstates of QCD in a finite cubic 
volume to the phase shift for elastic scattering. By ex- 
tracting multiple excited eigenstates within a given quan- 
tum number sector, one can map out the phase shift as a 
function of scattering momentum and, if present in that 
channel, observe resonant behaviour. 

In this letter, we demonstrate the feasibility of the 
technique in a simple sector, that of tttt scattering in 
isospin-2 (/ = 2), where the interaction is not strong 
enough to form a resonance, but rather is weak and re- 
pulsive. For the first time using this method, we extract 
the S- and D-w&ve phase shifts as a function of scatter- 
ing momentum. This procedure is carried out indepen- 
dently on multiple volumes to validate the finite-volume 
method. We find through computations at a range of 
quark masses that at the level of precision attained the 
phase shift is largely quark mass independent. 

Previous lattice QCD calculations of 7T7t scattering 
have limited themselves to extracting the phase shift at 
near-zero energy, more conveniently expressed via the 
scattering length [H [5] , or by performing the same cal- 
culation in a moving frame, a single phase shift point at 
non-zero energy [6]. 

In contrast, we use the "distillation" method [7] to con- 
struct both creation and annihilation operators of defi- 
nite 7T7T relative momentum, and employ them to form 
a variational basis of composite QCD operators that re- 
semble pairs of pions. This enables us to extract a spec- 
trum of multiple states with / = 2, £ p = + , 2 + (£ is the 
partial wave and P the parity) and, using the Liischer 
technique, we find the phase shift as a discrete function 



of the scattering momentum. This sets the groundwork 
for investigating resonances in other meson-meson scat- 
tering channels. 

Experimentally, 7T7t 1 — 2 phase shifts have been 
extracted from irN — > ttttN' charge-exchange scattering 
reactions, treating the dominantly-exchangcd pion as 
approximately on-shell owing to the proximity of the 
t-channel pole to the physical small-i region. The extant 
data [8HTT] for 5g= and Sg =2 are broadly consistent 
in the low-energy region measured and there is little 
statistically significant evidence for inelasticity. 

Finite volume analysis: Liischer's method relates the 
discrete spectrum of energy levels in a finite volume to 
phase shifts evaluated at the scattering momenta corre- 
sponding to the extracted energy values. Complications 
arise from the cubic symmetry of the lattice boundary 
which reduces the irreducible symmetry channels from 
the set of all integer spins to a finite set of irreducible rep- 
resentations. The relevant irreps, T, for tttt isospin-2 scat- 
tering at low momentum are which contains contin- 
uum spins £ = 0,4..., T+{1 = 2,4...), E+{£ = 2,4...) 
and {£ = 4 . . .) . Odd £ do not contribute due to Bose 
symmetry. 

Once the finite volume energy levels, E^ are ob- 
tained from an explicit Monte-Carlo calculation on a 
fixed volume (L 3 ) lattice, the scattering momenta fol- 
low assuming a continuum-like dispersion relation, k — 
\J (E 7T7T /2) 2 — to 2 . The desired phase-shifts are embed- 
ded in an equation 



det e 2i5 ^ 



0. 
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where Ur (k^) is a matrix in the space of partial waves, 
£, of known functions particular to this irrep, T, evalu- 
ated at the scattering momentum, k. e 2lS ^ is a diagonal 
matrix featuring phase-shifts, 5g{k), for all partial waves 
contributing to the irrep T. The dimension of these ma- 
trices is formally infinite, since there are an infinite num- 
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ber of possible partial waves contributing to each irrep 
r. However, one can argue that, since higher waves typ- 
ically contribute less at low momentum, one can cut-off 
the dimension at some low-spin and, provided the results 
are reasonably insensitive to this cutoff, reliably extract 
the phase shift for low partial waves. This is the core of 
the Liischer method [3]. 

The aim then is to solve for some set of phase-shifts 
{5i{k)}, but since Equation[l]is one equation (per energy 
level) in several unknowns, this will not be possible. 
Instead we will attempt to bound the size of all Sg(k) 
for I higher than the lowest in each irrep using other 
information. In practice we will assume that, in the 
energy region accessible to us, <5^>4 ~ and consider 
only the effect of a non-zero 64(h). 

Correlator construction and variational analysis: In 
order to obtain the finite-volume energy spectrum, we 
form a matrix of correlators using a basis of tttt operators 
that is then diagonalised using the variational method 
EH- The operators are constructed to resemble 
a pair of pions with total momentum zero and definite 
relative momentum: 

m p 

The subduction coefficients, <S r '™, project operators of 
definite I into definite irreps, T - their explicit forms can 
be found in Appendix A of The sum over directions 
of momentum, p, at a fixed magnitude is limited to those 
allowed by the periodic cubic boundary conditions. On 
a lattice with spatial extent L these are p = ^-n for a 
vector of integers n. 

In this first study we utilise only a simple operator 
capable of interpolating a pion at momentum p from the 
vacuum, 

oAp) = [^□,7 5 n^] (*), 

X 

where the quark fields are acted upon by a distillation 
smearing operator that emphasises the low momentum 
quark and gluon modes that dominate low mass hadrons. 
In this study we use D a = X) n =T e " l £n£j where 
A n , Cn are the eigenvalues and eigenvectors of the gauge- 
covariant three-dimensional Laplacian operator (see [21 E] 
for details; a = was used in [U [2] ) . It is distillation that 
factorises the construction of correlators in such a way as 
to make possible the projection onto definite inter-pion 
momentum at both source and sink, something that is 
not possible in the traditional "point-all" method. De- 
tails of the distillation correlator construction can be 
found in jJJ. 

Our variational basis in the irrep consists of opera- 
tors with \p\ 2 = (^-) 2 (0, 1, ... 4) each with two smearing 
radii a — 0.0, 4.0, giving a ten dimensional basis. For E + 
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TABLE I: Lattices used in this study. N VCCB indicates the 
number of eigenvectors of the laplacian used in the distillation 
method. 

we have \p\ 2 = (^-) 2 (1,2,4) and two smearings and T 2 + 

with |p| 2 = (77) 2 (2, 3) and two smearings. The T± irrep 
has lowest spin £ = 4 for two pions. However the lowest 
momentum from which a operator can be constructed 
is \p\ 2 = 5(^) 2 and this is the only one we used. With 
these operators at source and sink, we form all correlators 
using Wick contractions relevant for I = 2. 

Computations are performed on anisotropic lattices 
with three dynamical flavors of Clover fermions with spa- 
tial lattice spacing a s ~ 0.12 fm and finer temporal spac- 
ing, a^ 1 ~ 5.6 GeV, see Table [i] and [H]. A precise mea- 
sure of the anisotropy, £ = a s /at, is required to determine 
the spatial length of the lattice in temporal lattice units, 
L/at — ^L/a s . Fitting single-pion correlators at finite 
momentum, a s p — 2 J a n, determines atE v (\n\) and £ 
follows from fitting the dispersion relation 

{a t E^\n\)) 2 = (a t m^ 2 + ± \n\ 2 , 

for multiple values of \n\ and L/a s . Explicitly we 
find f = 3.459(4), 3.454(5), 3.459(3) on respectively the 
m v = 396,444,524 MeV lattices, showing the lack of 
quark mass dependence observed previously and utilised 
in the dynamical tuning of the lattice action[T5]. Mass- 
dimension quantities multiplied by the temporal lattice 
spacing, at, are scale-set using the procedure outlined in 
[2] , using the f2-baryon mass determined on the same lat- 
tice, m — atm ■ mo hys '. The continuum scaling of the 
results is not investigated in this calculation at a single 
lattice spacing. 

In Figure [l] we show the finite-volume spectra ob- 
tained with m w = 396 MeV and L/a s = 16,20,24. 
We clearly observe shifts relative to the energy of two 
non-interacting pions with back-to-back momentum of 
p = ^f-n, E-k-k = 2yjm 2 + \p\ 2 . It is this energy shift 
that Luscher's method relates to the scattering phase 
shift through Equation [l] It is the fact that we are able 
to resolve excited energy levels with a statistical preci- 
sion below 1% that makes possible an extraction of the 
scattering phase shift as a function of scattering momen- 
tum. 
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FIG. 1: Low-lying spectrum, in units of the temporal lat- 
tice spacing, of finite volume states in irreps Af,E + ,T^~ at 
= 396 MeV for L/a s = 16,20,24. The box height indi- 
cates the statistical uncertainty on the energy. Orange boxes 
correspond to states suspected of being nir* scattering states. 
Dashed lines indicate the non-interacting energy of pion pairs 
with the allowed lattice momenta between them. 



Phase shift: For each tttt energy level in each irrep 
in Figure [l] we can set up an Equation [l] to be solved 
for the phase shifts, Sg. The simplest way to solve these 
equations is to neglect the contribution of £ > 4 to obtain 
<!>o from and 82 from or E + . Doing so gives the 
red, green and blue colored points in Figure [2j The small 
discrepancies between E + and extractions of 82 at 
k 2 ~ 0.35, 0.55, 0.85 GeV 2 (corresponding to the levels at 
a t E ~ 0.25, 0.29, 0.36 in Figure [l]), have a possible origin 
in the neglect of a non- negligible value of 84. We can 
estimate the size of this 84 by solving the coupled system 
of Equations [I] for T% and E + at the relevant energy for 
the two unknowns, 82,84. The values of 84 so extracted 
are shown by the pink points in Figure [2] 

For a direct estimate of 84 from T^, only the L/a s — 24 
lattice has a point within our plotted range of scattering 
momentum. The extracted point is shown by the pink 
diamond in Figure 2, and is in good agreement with the 
other estimates, showing that |#4 1 is less than 2° over the 
whole of the explored momentum range. 

With an estimated magnitude of 84(h) in hand (from 
interpolation between the determined points), we can 
solve Equation [T] including the effect of the £ — 4 wave. 
This gives rise to the orange, light green and cyan colored 
points in Figure [2] which are seen to differ relatively little 
from the points with ^4 assumed to be zero. For final 
presentation we enlarge the errorbar to include the effect 
of the estimated 84 giving rise to asymmetric errorbars 
in Figure [3j 

As indicated in Figure [2] the 47r threshold opens 
within the energy range of our extracted phase-shifts 
and technically for energies above this the formalism 
leading to Equation [l] is not rigorously correct. On the 
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FIG. 2: Phase-shifts extracted from spectra with = 
396MeV. Red (A+), green (E + ), blue (T 2 + ) colored points 
assume £4 = 0; orange (A^ ), light green (E + ), cyan (T^~) 
colored points (shifted slightly to the right) used estimated 
£4 as described in the text: note that the corrected 82 val- 
ues from E + ,T^ coincide by construction at momenta near 
\p\ 2 = 2 ■ f^M • Estimated $4 shown by pink points. Also 
indicated are the positions of inelastic thresholds into 4vt and 
pp. 



other hand, there is relatively little evidence experimen- 
tally for considerable inelasticity in the 7T7t isospin-2 
channel in the energy range so-far probed - what little 
data there is does not show statistically significant devi- 
ation from an elastic approximation \Q . As an initial 
approximation, we shall assume that the inelasticity is 
negligible and continue to use Equation [T] above the 
inelastic threshold. The £ = 2 phase-shift extracted from 
E + , T2 irreps should be less sensitive to any inelasticity 
since the effective threshold in finite- volume is higher as 
it requires at least one unit of relative momentum in the 
47T system. Future calculations should test the elasticity 
assumption by computing correlators using operators 
that resemble four pions projected into isospin-2 in the 
appropriate partial waves. 

Results: In Figure [3] we show our results for S and D- 
wave phase shifts, at a range of pion masses, along with 
experimental data taken from [SMllj. We observe rea- 
sonable agreement with the experimental data at lower 
scattering momenta, where the scattering is purely elas- 
tic, for all the pion masses computed. This suggests that 
it is possible that the phase-shift is only mildly dependent 
upon pion mass. Of course, one requires lattice computa- 
tions at smaller pion masses to verify that the agreement 
with experiment continues. 

Using only 7T7t and 7r correlators at zero momentum 
we can perform the extraction of the scattering length 
using the methodology of [H |3] . The scattering lengths 
so obtained are shown in Figure |4j where they are seen to 
be in reasonable agreement with the precision data of [1] , 
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FIG. 3: Phase-shift in degrees for tttv 1 = 2 scattering with 
I = 0(<5o) and £ = 2(52). Lattice results at various pion masses 
and volumes. Experimental data from [8Tlll|. 



Summary and prospects: We have demonstrated the 
feasibility of an explicit application of the Liischer finitc- 
volume framework in dynamical lattice QCD. Using mul- 
tiple excited state energy levels extracted in a single vol- 
ume, we have determined the S- and D-wave tttt isospin-2 
phase shifts as a function of scattering momentum. Mul- 
tiple volumes are then used for validation and estimation 
of the effect of neglected higher partial waves. We esti- 
mate that |<5 4 | < 2° for k < 1 GeV. 

We observe no significant pion mass dependence in the 
phase shift below k ~ 1 GeV, with results for m, > 
400 MeV being in reasonable agreement with experimen- 
tal data at low scattering momentum. For precision cov- 
erage of experimentally relevant kinematics, we would 
require still larger volumes to sample points at smaller 
scattering momentum in the elastic region. 

This calculation sets the groundwork for an investiga- 
tion of the resonances in meson-meson scattering that 
arise from the strong interaction. Inclusion of quark 
annihilation diagrams in the calculation of correlators 
[17] will enable the I = 1 tttt sector to be studied, where 
one expects to see the p resonance appearing as a rapidly 
rising phase shift. Some attempts in this direction have 
been made [18 - 20 , but using only a small basis of 
operators and subsequently extracting a very limited 
number of phase-shift points. Distillation and stochastic 
variants [211 E2] will allow us to efficiently construct 
a large basis and thus map out many points on the 
phase-shift curve. In future work, we will explore a range 
of different scattering hadrons in various partial-waves. 



computed on a lattice of similar spatial lattice spacing. 

We can also obtain estimates for the scattering length 
and effective range by fitting the k dependence of So(k), 
where we find that scattering lengths largely agree with 
the estimates from the simple method above, while the 
effective range is small but only poorly determined. 
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FIG. 4: 5-wave scattering length. Blue points from p = 
correlators, pink points from effective range fits to S(k). Com- 
parison to lattice results of [U [5] and Roy equations analysis 
of experimental data [16] . 
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